Next-to-leading order corrections for gg ? ZH with top quark mass dependence

In this Letter, we present for the first time a calculation of complete next-to-leading order corrections to $gg \to ZH$ process. We use method small mass expansion tackle most challenging two-loop virtual amplitude, in which top quark dependence is retained throughout calculations. show that our provides reliable numeric results all kinematic regions, and phenomenological predictions total differential cross sections at Large Hadron Collider its future upgrades. Our are necessary ingredients towards reducing theoretical uncertainties $pp down percent-level, provide important inputs precision experimental collider programs.